Pairing and the structure of the p/-shell N ~ Z nuclei 
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The influence of the isoscalar and isovector L = pairing components of the effective nucleon- 
nucleon interaction is evaluated for several isobaric chains, in the framework of full pf shell model 
calculations. We show that the combined effect of both isospin channels of the pairing force is re- 
sponsible for the appearance of T = 1 ground states in N — Z odd-odd nuclei. However, no evidence 
is found relating them to the Wigner energy. We study the dependence of their contributions to the 
total energy on the rotational frecuency in the deformed nucleus 48 Cr. Both decrease with increasing 
angular momentum and go to zero at the band termination. Below the backbending their net effect 
is a reduction of the moment of inertia, more than half of which comes from the proton-neutron 
channel. 



The study of the isovector pairing interaction among 
like particles is one of the classical themes of nuclear 
physics. Proton-neutron pairing, has been much less 
studied, in particular its isoscalar part. Early attempts 
to solve the isoscalar and isovector L = pairing in one I 
shell can be found in ref. Q and an extension to the two 
shells case in [^). The experimental access to medium 
mass N = Z nuclei, approaching the proton drip line has 
renewed the interest for the components of the nuclear 
interaction whose manifestations are enhanced close to 
N = Z. The proton-neutron part of the isovector pair- 
ing near N = Z has been studied in term of algebraic 
model and compared with Shell Model Monte Carlo cal- 
culations in ref. S. 

In particular there are many discussions on the role of 
the isoscalar pairing interaction. Is it responsible for the 
Wigner energy? Should it give rise to a proton-neutron 
condensate? How will its manifestations evolve as the an- 
gular momentum increases?. Some of these issues have 
been examined in refs. The role of the isoscalar 

pairing channel in the processes mediated by the Gamow- 
Teller operator has been recently studied in ref. ||. 

In this letter we take advantage of the availability of 
detailed microscopic descriptions of the nuclei in the mid- 
dle of the p/-shell -that have already passed many ex- 
perimental tests- to address some of these questions j^U • 
The existence of well deformed rotors in the A = 50 mass 
region makes it possible to explore the dependence of a 
given observable with the rotational frequency || |T^| . 

We shall study the isobaric chains A — 46, 48 and 50, 
using the spherical shell model approach. The details of 
the shell model work have been explained in ref. (?j . We 



use the Lanczos shell model code ANTOINE (ref. (13 
and the the effective interaction KB3. All the calcula- 
tions presented are unrestricted in the p/-shell. 

A key ingredient of our approach is the choice of the 
"physical" isocalar and isovector pairing terms of the ef- 



fective nuclear interaction. We follow here the work of 
Dufour and Zuker in ref [Jh|, where they give a precise 
characterization of the multipole part of the nuclear in- 
teraction, showing in particular that the pairing term of 
any realistic G-matrix has only two contributions; the 
well known schematic isovector L = 0, S = pairing 
and an isoscalar part which turns out to be again the 
schematic L = 0, S = 1 pairing. This eliminates many 
ambiguities associated to the choice of the isovector and 
isoscalar pairing interactions. We shall keep the nota- 
tion in E3] and use P01 for the isovector and P10 for the 
isoscalar pairing. The explicit expressions for the two 
body matrix elements in LS and jj couplings are: 

(I a IbLSJT\P01\I c IdLSJT) 

= <V(2U + l)(2Zc + 1)SabScdS L0 S S0 Sj Sti, (la) 
(I a IbLSJT\PW\I c IdLSJT) 

= G^(2l A + 1)(2/d + 1)5 A bScdSloSsiSji6to, (lb) 



(jaj b JT\P01\ JcJd JT) 

= G^(j a + 1/2) (j c + l/2)6 ab 5 cd 8j 8 T i, 

(j a jbJT\P10\j c j d JT) 
2(_i)jWc 



(lc) 



G 



V(2ja + l)(2j 6 + l)(2j c + l)(2j d + 1) 



1/2 j a l a 
Jt 1/2 1 



1/2 Jc lc 
3d 1/2 1 



SabScdSjxSto, 



(Id) 



where capital letters denote a /-orbit (quantum numbers 
nl) and lowercase letters denote a j-orbit (quantum num- 
bers nlj). The value for the pairing constant G is easily 
obtained from the numbers given in table I of the same 
reference. We use G = -0.295 for P01 and G = -0.459 
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for P10. It could be thinkable to do the same kind of 
analysis with a monopole free pairing interaction. We 
have verified that none of the conclusions of the paper 
would be modified if we did so. Thus, we decided to 
stick to the usual forms of pairing. 

Once this definition of the pairing operators adopted, 
we focuss in the role they play in the behaviour of differ- 
ent nuclei and different physical quantities. For that we 
shall proceed in the most straightforward way; for each 
case we make first a reference calculation using the inter- 
action KB3, then we substract from KB3 the isovector or 
the isoscalar pairing hamiltonians and make the calcula- 
tions with the new interactions KB3-P01 and KB3-P10. 
We obtain the effect of each pairing channel by direct 
comparison with the reference calculation. One could 
also estimate these effects in perturbation theory, in the 
cases studied here we have verified that the results are 
equivalent, however this may not be true in general. 



Binding energies for an isobaric chain. 

We begin with A — 46 in figure |[ We have plot- 
ted the contributions to the ground state binding energy 
(BE) from the isovector pairing (labeled POI) and from 
the isoscalar pairing (labeled P10). BE's are taken with 
the positive sign. The figure shows the strong odd-even 
staggering of the POI points, mostly suppressed in the 
P10 expected. The effect of PIO, although im- 

portant, is smaller than that of POI and goes smoothly 
to zero as the number of valence protons goes to zero 
(notice that our description is fully simetric under the 
exchange of protons and neutrons). The only little sur- 
prise is that moving from T = 1 to T = 0, not only the 
contribution of POI decreases, but also the contibution of 
P10, although to a lesser extent. This may explain why 
the ground state of 46 V has T = 1 instead of the expected 
T = 0. It goes like this; the monopole part of KB3 will 
put the centroid of the T = states lower than the cen- 
troid of the T — 1 states by about 1.3 MeV; on the other 
hand, the total pairing contribution to T — 1 is larger by 
more than 1.5 MeV than the contribution to T = (see 
figure [j]), therefore it is finally the T = 1 state that be- 
comes the ground state of the odd-odd T z = member of 
the multiplet. Notice that the inversion depends on the 
balance between the isovector monopole term and both 
pairings. As the isovector monopole has a smooth de- 
pendence with mass, while the isovector pairing depends 
linearly with the degeneracy of the orbit being filled, we 
expect the inversion to take place when high j orbits 
are dominant. When low j orbits are the relevant ones, 
the isoscalar pairing contribution to the T — state can 
be enhanced due to the small spin-orbit splitting, what 
helps producing T = ground states. This is borne out 
by the experiments; when the I/7/2 orbit is being filled 
the T = 1 states are ground states, in 58 Cu the 2p 3 / 2 
orbit is dominant and restores a T = ground state, in 
62 Ga the influence of the I/5/2 orbit shifts again to a 



T = 1 ground state. From there on the I59/2 orbit starts 
to play an important role and T — 1 ground states occur 
in all cases. Actually the T — 0-T = 1 splitting could 
give a hint on the ordering of the single particle orbits in 
places where no other information is available. 

In figure || the same quantities are plotted for A = 48. 
The behaviour of the P10 contribution is very smooth 
and we do not find any breaking of slope at T = 0. 
The POI points stagger as usual. In figure || we present 
A = 50. Now the part of the isobaric chain lying in- 
side the p/-shell is longer and the trends seen in A = 46 
better established. There is a slight increase of the P10 
contribution from T — to T = 1 followed by a smooth 
decrease, reaching zero at T — 5. As in A = 46 the re- 
duction of both P10 and POI contributions from T = 1 
to T = explains the occurrence of a T = 1 ground state 
in 50 Mn. 



Wigner energy. 

This name is given to the extra binding of the T = 
member of an isobaric multiplet compared to what would 
be predicted by an (N — Z) 2 extrapolation from the other 
isobars. We can make an assesment of the role of the POI 
and P10 terms proceeding as in the definition just given. 
In A = 46 we have too few nuclei in the valence space as 
to make the fit, therefore we only present results for the 
two other chains. 
For A = 48, 

e w = 2.82 KB3 

e w = 3.04 KB3-P01 

e w = 2.47 KB3-P10 

clearly the effect of both types of pairing on the Wigner 
energy is negligible. 
For A = 50, 

e w = 1.99 KB3 

e w = 2.60 KB3-P01 

e w = 2.12 KB3-P10 

the effect of the pairing terms is even less significant than 
in the A = 48 isobars. In conclusion, we find no link be- 
tween the Wigner energy and the dominant pairing terms 
of the nuclear interaction. Furthermore, if we had made 
the fit with the more microscopic prescription T(T + 1) 
the values of this "new" Wigner energy would have been 
much smaller (1.23 and 0.24 for A = 48 and A = 50) and 
the effect of the pairing terms negligible again. We there- 
fore think that the Wigner energy cannot be explained 
as a pairing effect. 



Pairing and angular momentum 

The evolution of the pair content of a nucleus as the ro- 
tational frequency increases has been the subject of many 
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heated debates in the recent past. We have studied this 
evolution in several cases but we shall present only re- 
sults for 48 Cr, which is the most representative example. 
In figure |] we plot the contributions of the P10 and POI 
pairing hamiltonians to the absolute energy of the yrast 
levels. Two regions show up, corresponding roughly to 
below/above the backbending. In the region of collective 
rotation the POI contribution, that at J = almost dou- 
bles the P10, decreases with J much more rapidly than 
the P10 one. At J = 8 both reach approximately the 
same value, and stay constant until they suddenly vanish 
at the aligned J = 16 state. In figure || we have trans- 
lated the results of figure |4| to a backbending plot. The 
most visible effect of the pairing correlation is the well 
known reduction of the moment of inertia. The effect of 
the POI channel is clearly dominant but the P10 contribu- 
tion is far from being negligible and amounts to about one 
fourth of the POI. This suggest a plausible explanation to 
the discrepancies found in the comparison of the Cranked 
Hartree Fock Bogolyubov results using the Gogny force 
and the Shell Model (or the experimental ones) for the 
yrst band of 48 Cr ||. Both approaches were shown to 
predict correctly the quadrupole properties, however the 
CHFB calculation gave a much larger moment of inertia. 
Actually it came out a factor three larger than the exper- 
iment, which is fully compatible with what we have seen 
to be the effect of taking out the proton neutron pair- 
ing (T = and T = 1). Indeed, proton-neutron pairing 
is not properly taken into account in the standard HFB 
treatements. 50 Cr and 50 Mn behave in a qualitatively 
similar way. The only difference being the smaller role 
played by the P10 pairing in the moment of inertia of 
50 Cr. 



Pairing condensates? 

There has been speculations recently on the possible 
existence of a new superfluid phase of nuclear matter 
made of a condensate of L = 0, S = 1 neutron proton 
pairs, that could be approachable in N = Z nuclei. Once 
the meaning of "condensate" agreed, we can try to con- 
tribute to the debate. 

Imagine that the I/7/2 orbit were the valence space 
and POI the effective interaction, in this case, the ground 
state of an even-even nucleus will be the seniority zero 
state that can be easily interpreted as a condensate of 
L = 0, 5 = 0, T = 1 pairs (with all the caveats, because 
only the valence particles are paired and the j's of the 
orbits are not large enough compared with the number of 
active particles). Given that the (1/7/2)™ configurations 
are the leading ones in the ground states of the nuclei 
we are interested in, it is in principle possible that the 
"POI condensate", i.e. the ground state of the (1/7/2)™ 
isovector pairing problem, and the physical ground state 
could have a reasonable overlap. Furthermore, most of 
the properties of the (I/7/2)™ pairing condensate persist 
if we solve the pairing problem in the full p/-shell with 



the experimental single particle splitting. 

On the contrary, if we solve the hamiltonian P10 in 
the (I/7/2)" configurations, we do not find any conden- 
sate. In order to obtain a condensate of L = 0, S = 1, 
T = pairs, the valence space must include at least the 
I/5/2 orbit degenerate with the I/7/2, or, in general to be 
made of Z-orbits, which is an extremely unphysical situa- 
tion. Taking this argument to its limit one could say that 
the existence of a triplet pairing condensate is excluded 
by the spin orbit splitting of the nuclear mean field. 

We have tried to make these statements quantitative 
computing the overlaps of the wave fuctions of the two 
condensates and the physical ground state in several pf- 
shell nuclei. Let's call |KB3) the physical ground state 
of a given nucleus obtained with the interaction KB3; 
|CP01) the isovector condensate obtained solving POI in 
the pf shell with the experimental single particle split- 
ting and |CP10) the isoscalar condensate obtained solv- 
ing P10 in the pf shell without spin-orbit splitting, i.e. 
I = 3 and I = 1 separated by 2.0 MeV. The squared over- 
laps of these states with the physical ground state give 
the probabilities of finding the nuclei in the different con- 
densed phases. The results are gathered in table I and the 
conclusions rather evident. The isovector L — pairing 
condensate clearly dominates the structure of the Cal- 
cium isotopes we have selected. When two protons are 
added the overlaps get reduced and the interpretation in 
term of condensates is less justified. For the Chromiums 
no trace of a condensate is left. Notice that an overlap 
of 38% can be considered almost marginal. Actually, the 
overlap of the 48 Cr ground state, obtained using a pure 
Q ■ Q interaction and the isovector pairing condensate 
is already 25%. As anticipated the L=0 isoscalar pair- 
ing condensate does not show up, on the contrary its is 
almost exactly orthogonal to the physical ground states 
in most of the cases. As we said before, the existence of 
such a condensate would only be possible if in the valence 
space that accomodates the physics of the problem, the 
spin-orbit splittings were negligible. 



Conclusions 

We have explored the effect of the isovector and 
isoscalar L — pairing interactions on several properties 
of p/-shell nuclei close to N — Z. We have found that the 
isoscalar pairing, whose contribution to the ground state 
binding energy of the nuclei studied may be very impor- 
tant, is not at the origin of the Wigner energy. Nor the 
isovector term. We argue that the occurrence of T = 1 
ground states in odd-odd, N = Z nuclei is due to the 
combined action of the two pairing channels that, when 
the dominant orbit has large I, can overcome the isovector 
monopole gap. We stress the need of a good treatment of 
the proton-neutron pairing in order to obtain the correct 
moment of inertia oi N = Z rotors. Finally we discard 
the presence of isoscalar pairing condensates everywhere 
in the region, as well as the existence of isovector pairing 
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condensates at N = Z. 
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TABLE I. Overlaps between the physical ground state and 
the singlet and triplet pairing "condensates" in several nuclei. 

nucleus |(KB3|CP01)] 2 (KB3|CP10) \ 2 |(CP10|CP01)| 2 



Ca 97 - 

Ca 95 

; Ti 63 16 1 

; Ti 52 10 1 

; Ti 59 1 

Cr 38 2 

Gr 42 1 
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FIG. 1. Pairing contributions to the ground state energies 
of the A=46 isobaric multiplet (in MeV) 



FIG. 3. Pairing contributions to the ground state energies 
of the A=50 isobaric multiplet (in MeV) 
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FIG. 2. Pairing contributions to the ground state energies 
of the A=48 isobaric multiplet (in MeV) 



FIG. 4. Pairing contributions to the energies of the yrast 
states in 48 Cr. 
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FIG. 5. Gamma ray energies along the yrast band of 48 Cr 
(in MeV); full interaction (KB3); without isoscalar pairing 
(KB3-P10); without isovector pairing (KB3-P01) 
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